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SUMMARY 


By  the  use  of  a  form  of  the  Prandtl-Glauert  method  that  is  valid 
for  three-dimensional  flow  problems,  the  value  of  the  maximum  Incremental 
velocity  for  compressihle  flow  about  thin  ellipsoids  at  zero  angle  of 
attack  is  calculated  as  a  function  of  the  Mach  number  for  various  aspect 
ratios  and  thickness  ratios.  The  critical  Mach  numbers  (within  the 
accuracy  of  the  Prandtl-Glauert  method)  of  the  various  ellipsoids  are 
also  determined.  The  results  indicate  an  Increase  in  critical  Mach 
nimber  with  decrease  in  aspect  ratio  which  is  large  enough  to  explain 
experimental  results  on  low-aspect -ratio  wings  at  zero  lift. 


INTRODUCTION 


Recent  tests  (references  1  and  2)  have  shown  that  an  appreciable 
increase  in  the  critical  Mach  number,  together  with  other  improvements 
of  the  aerodynamic  characteristics  at  supercritical  Mach  numbers,  results 
from  the  use  of  wings  of  very  low  aspect  ratio.  These  Improved  charac¬ 
teristics  have  been  somewhat  q.ualitatively  ascribed  to  "three-dimensional 
relief,  although  no  q.uantltatlve  theoretical  discussion  has  yet  been 
provided. 

In  the  present  paper  an  effort  is  made  to  provide  such  a  study  by 
considering  the  flow,  at  zero  angle  of  attack,  about  a  series  of  thin 
ellipsoids  of  various  aspect  ratios  and  thickness  ratios.  Ellipsoids 
were  chosen  because  they  are  amenable  to  calculation.  Although  they 
differ  appreciably  from  the  wings  of  reference  1,  which  had  an  NACA 
0012  airfoil  section  and  rectangular  plan  form,  ellipsoids  should  never¬ 
theless  show  similar  aspect -ratio  effects.  The  calculations  were  made 
for  ellipsoids  of  thickness  ratios  0.10,  0.15^  and  0.20,  and  for  the 
entire  range  of  aspect  ratios  from  the  elliptic  cylinder  to  the 
ellipsoid  of  revolution. 

The  compressibility  effects  were  computed  by  the  use  of  a  form  of 
the  Prandtl-Glauert  method  that  is  valid  for  three-dimensional  flow 
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problems.  The  method  has  heen  given  hy  G8thert  (reference  3)  without, 
however,  very  explicit  mathematical  proof.  Another  correct  statement 
of  the  three-dimensional  form  of  the  Prandtl-Glauert  method  was  given 
earlier  by  A.  Busemann  In  reference  4,  where,  however,  no  formulas  were 
given.  Since  the  methods  that  have  been  commonly  used  (see,  for  example, 
references  5^  6,  and  7)  are  applicable  only  to  two-dimensional  problems, 
a  detailed  proof  of  the  method  correct  for  three-dimensional  flow  Is 
Included  In  the  appendix.  A  brief  discussion  of  the  accuracy  of  the 
Prandtl-Glauert  method,  as  applied  to  ellipsoids,  is  also  given.  This 
work  was  completed  in  April  1946. 


SYMBOLS 


U 

C 

M 

7 

3  =  \|l  -  m2 

B 

9 

u,  Y,  w 

B' 

U*,  V*,  w' 

a 

b 


free-stream  velocity 

velocity  of  sound  in  free  stream 

free-stream  Mach  nimiber  (u/c) 

ratio  of  specific  heats  (7  =  1.4  for  air) 


rectangular  coordinates 
thin  body 

velocity  potential 

X-,  y-,  and  z-components  of  incremental 
velocity  for  compressible  flow  about  B 

body  obtained  by  stretching  B  in  direction 
of  x-axls  by  the  factor  l/p 

X-,  y-,  and  z-components  of  incremental 

velocity  for  incompressible  flow  about  B ’ 

maximum  semi  chord  of  ellipsoid 

semi span  of  ellipsoid 


c 


maximum  semi  thickness  of  ellipsoid 
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A 


^  ”  U 

u(M) 

u(0) 


aspect  ratio 


value  of  u  vhen  the  Mach  number  is  eq.ual  to  M 
value  of  u  for  incompressible  flow  (M  =  O) 


€ 


thickness  ratio 


(Thickness 
Chord  ^ 


(€,M)  value  of  ratio  of  Incremental  velocity  to 

free-stream  velocity  for  compressible 
flow  having  Mach  number  M  about  a  body 
having  thickness  ratio  6 


(e,0)  value  of  ratio  of  incremental  velocity  to  free- 

stream  velocity  for  incompressible  flow  about 
a  body  having  thickness  ratio  € 

Subscript: 

max  maximum  value 


METHODS  OF  CALCULATION 


The  Frandtl-Glauert  method  for  three-dimensional  flow*  -  The 
Prandtl-Glauert  method  is  used  in  the  present  paper  in  the  following 
form: 


The  Incremental  velocities  at  a  point  P  of  a  three-dimensional 
compressible  flow  field  about  a  thin  body  B  may  be  obtained  in  three 
steps: 

(1)  The  x-coordlnates  of  all  points  of  B  are  Increased  by 
the  factor  l/p,  where 

3  =  \]l  -  m2 

and  where  the  x-axis  is  in  the  stream  direction.  This  transformation 
takes  B  into  a  stretched  body  B'. 
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(2)  The  incremental  velocities  u*,  v*,  w*,  in  the  direction  of 

the  j-y  and  z-axes^  respectively,  at  the  point  P*  in  the  flow 

field  of  B  *  corresponding  to  the  point  P  in  the  flow  field  of  B 
are  calculated  as  though  B*  were  in  an  incompressihle  flow  having  the 
same  free-stream  velocity  as  the  original  compressihle  flow- 

(3)  The  values  u,  v,  and  w  of  the  incremental  velocities  at  the 
point  P  in  the  compressihle  flow  field  of  the  original  imstretched 
hody  are  then  found  hy  the  ecLuations 


1  , 
w  =  —  w 

P 


A  derivation  of  this  form  of  the  Prandtl-Glauert  method  is  given 
in  the  appendix.  The  method  in  essentially  this  form  has  heen  given  hy 
Gothert  (reference  3)  without,  however,  a  very  clear  proof.  (Go'thert 
prefers  to  shrink  the  lateral  coordinates  of  the  hody  hy  the  factor  P 
rather  than  to  expand  the  coordinate  in  the  stream  direction  hy  the 
factor  l/P;  obviously  the  two  procedures  lead  to  the  same  result.) 

Prandtl  (reference  5)  find  Von  Karman  (reference  6)  state  the  method  in 
a  form  that  is  valid  for  two-dimensional  flows  hut  in  general  is  incorrect 
for  three-dimensional  flows.  Goldstein  and  Young  (reference  7)  filso  give 
a  discussion  leading  to  results  that  are  correct  only  for  two  dimensions. 

A  discussion  of  the  reasons  for  the  failure  of  these  commonly  used 
methods  for  three-dimensional  flow  prohlems  is  included  in  the  appendix. 

Calculation  of  incremental  velocity  for  comnressihle  flow  about 
ellipsoids.  -  In  order  to  determine,  hy  the  Prandtl-Glauert  method,  the 
incremental  velocity  on  the  surface  of  an  ellipsoid  having  semiaxes  a, 
h,  and  c.  (where  a  ■  is  the  length  of  the  semi  axis  in  the  stream  direc¬ 
tion)  ,  the  Incremental  velocity  is  calculated  for  a  stretched  ellipsoid 

having  semiaxes  a*,  h,  and  c  (where  a'  =  ^  in  an  incompressihle 
flow  having  the  same  stream  velocity)  and  the  result  is  multiplied 

hy  1/3^.  For  incompressible  flow  about  the  stretched  ellipsoid,  the 
velocity  potential  on  the  surface  of  the  ellipsoid  is  given  hy 


a. 


a. 


Ux 
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•where 


a. 


_ dX  _ 

(^a'2  +  x}  ^(a'2  +  x)  (y^  +  X^  Qi^  +  x) 


(See,  for  example,  reference  8,)  The  Incremental  velocity  at  x  =  0 
(half-chord  line  on  the  stretched  ellipsoid  In  Incompresslhle  flow)  Is 
then  given  hy 


u'  = 


U 


This  value  Is  the  maximum  value  of  u’  (reference  8)  and  evidently 
Is  the  same  at  all  points  on  the  half -chord  line*  The  Incremental 
velocity  at  the  half-chord  line  for  the  compressible  flow  about  the 
original  unstretched  ellipsoid  Is  given  by 


u 


1 

3^  2  -  tto 


U 


(1) 


Various  formulas  are  necessary  for  the  evaluation  of  the  Integral 
when  a*>b>c,  b>a^>c,  or  a'>b  =  c  (ellipsoid  of  revolution). 

For  a'  >  b  >  c,  the  value  of  Is  given  by  the  formula 

“o  -  ; - —  , .  (F  -  K)  (2) 

(a'2  -  1,2) 

where  F  and  E  are  incomplete  elliptic  integrals  of  the  first  and 
second  kind,  respectively,  defined  as  follows: 


Jo  Vl  -  sln^ 


_ _ 

yl  ■  slii2i|f  d-^f 


E  = 
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A 


where 


and 


»2  ~ 

sin  <p  =  - j - 


For  h  >  a*  >  c,  the  value  of  a  is  given  hy  the  formula 


ga'bc 

-  a  •2)  (a'=  -  o") 

E  are  defined  hy  the  foregoing  formulas  hut 

-  a '2 

Vb2  -  c" 


h 

Eq.uation  (2)  is  derived  from  the  first  equation  given  in  equa¬ 
tions  (5*13)  of  reference  8  hy  substituting  a*  for  a  and  hy  using 
the  expression  for  k  in  terms  of  a',  h,  and  c-  Equation  (3)  is 
derived  from  the  second  equation  given  in  equations  (5*13)  of  reference 
hy  interchanging  a  and  h,  substituting  a’  for  a,  and  using  the 
expression  for  k  in  terms  of  a*,  h,  and  c. 

For  a’  >  h  =  c  (ellipsoid  of  revolution),  aQ  is  given  hy  the 
equation 


and 


k  = 


sin  9  = 


G 


where  F  and 


ttg  =  a 


_ dX _ 

(a '2  +  x2)3/2(b2  +  x) 
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which  resolves  Into 


where 


If  this  value  for  cLq  is  substituted  in  ec[uation  (l) ,  the  incremental 
velocity  at  the  half-chord  line  for  the  ellipsoid  of  revolution  is  found 
to  be  . 


r 


u  = 


1 

32 


log, 


1  +  e 
e  1  -  e 


-  2e 


2e 


1  -  e' 


“  log. 


1  +  e 
1  -  e 


U 


The  limiting  case  of  infinite  aspect  ratio  (elliptic  cylinder)  was 
treated  by  the  use  of  formulas  for  the  ellipse  in  two-dimensional  flow 
(reference  9)* 

Calculation  of  the  critical  Mach  number. -  For  flow  about  a  two- 
dimensional  body,  the  free-stream  Mach  number  for  which  sonic  speed  is 
first  reached  at  some  point  on  the  surface  is  called  the  critical  Mach 
number,  because  of  the  development  of  a  shock  and  the  accompanying 
deterioration  of  the  aerodynamic  characteristics  shortly  after  this 
Mach  number  is  exceeded.  If  an  extension  of  the  definition  of  the 
critical  Mach  number  to  the  general  three-dimensional  body  is  desired, 
the  definition  appears,  at  first  sight,  to  req.uire  some  revision,  since 
for  the  general  three-dimensional  case  the  shock  formation  on  parts  of 
the  body  may  occur  along  lines  yawed  with  respect  to  the  free-stream 
velocity  (reference  lO) .  The  boundary  lines  of  the  supersonic  regions 
(sonic  lines)  must,  however,  always  contain  a  portion  that  is  normal 
to  the  free-stream  velocity  and  thus  the  definition  of  the  critical 
Mach  number  for  the  two-dimensional  case  (infinite  unyawed  cylinder) 
may  be  extended  without  modification  to  the  general  three-dimensional 
case.  For  the  special  three-dimensional  case  of  the  infinite  yawed 
cylinder,  the  portions  of  the  boundary  lines  of  the  supersonic  regions 
that  are  normal  to  the  free-stream  velocity  are  represented  by  two 
points  at  infinity,  only  one  of  which  (at  the  downstream  end)  has  the 
necessary  qualities  for  accumulating  disturbances,  that  is,  for  shock 
formation.  It  can  be  seen  that  for  the  general  three-dimensional  body 
the  critical  Mach  number,  ^though  still  defined  in  the  same  manner  as 
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for  th^  two-dimensional  tody,  largely  loses  its  critical  significance, 
since  the  shock  may  "begin  to  form  over  only  a  very  small  part  of  the 
surface,  so  that  its  occurrence  does  not  necessarily  imply  an  imminent 
deterioration  of  the  aerodynamic  characteristics  of  the  "body* 

For  the  special  case  of  the  unyawed  ellipsoid  considered  in  the 
present  paper,  however,  no  appreciable  analysis  of  shock  foimation  or 
shock  extent  along  the  lines  Just  Indicated  seems  to  he  req_uired«  As 
is  shown  in  the  section  "Calculation  of  the  Incremental  Velocity  for 
Compressible  Flow  about  Ellipsoids",  the  maximum  velocity  for  an  unyawed 
ellipsoid  is  in  the  stream  direction  and  occurs  simultaneously  at  all 
points  along  the  half -chord  line.  Sonic  velocity  is  thus  reached 
simultaneously  along  a  line  that  extends  across  the  entire  span  of  the 
body  and  is  normal  to  the  stream  direction.  These  conditions  also 
exist  in  the  case  of  the  unyawed  infinite  cylinder,  that  is,  the  two- 
dimensional  body. 

The  critical  Mach  number  of  the  ellipsoid,  within  the  accuracy  of 
the  Prandtl-Glauert  method,  was  accordingly  determined  by  solving 
graphically  the  equation 


where  u(m)  Is  the  ratio  of  the  Incremental  velocity  at  the  half -chord 
line  to  the  stream  velocity  at  the  Mach  number  M- 

Accuracy  of  the  Prandtl-Glauert  method.  -  The  Prandtl-Glauert  method 
is  based  on  the  assumption  of  small  perturbations.  Consequently,  near 
the  nose  of  the  ellipsoids  discussed  in  the  present  paper,  where  the 
assumption  of  small  perturbations  is  violated,  the  results  given  by  the 
Prandtl-Glauert  method  cannot  be  expected  to  "be  reliable.  More  reliable 
values,  however,  should  be  obtained  for  the  maximum  Incremental  velocity, 
which  occurs  at  the  half -chord  line.  The  accuracy  of  the  Prandtl-Glauert 
approximation  for  the  maximum  incremental  velocity  may  be  estimated  by 
comparison  with  more  exact  solutions  of  the  compressible  flow  problem. 

An  iteration  method  in  which  the  Prandtl-Glauert  method  is  used  as  the 
first  approximation  has  been  proposed  by  Busemann  (reference  ll) .  The 
first  and  second  approximations  have  been  calculated  by  Hantzsche  and 
Wendt  for  the  elliptic  cylinder  (reference  12)  and  by  Schmieden  and 
Kawalki  for  the  ellipsoid  of  revolution  (reference  13)*  Calculation 
of  the  maximum  Incremental  velocity  for  the  elliptic  cylinder  having 
thickness  ratio  0.20  by  a  formula  for  the  second  approximation  given 
in  reference  11  shows  that  the  value  given  by  the  Prandtl-Glauert  method 
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at  a  Mach  number  of  0-8  is  almost  20  percent  lower  than  the  value  given 
by  the  second  approximation.  For  the  ellipsoid  of  revolution,  however, 
the  value  of  the  maximum  Incremental  velocity  given  by  the  Prandtl-Glauert 
method  agreed  with  the  value  given  by  the  second  approximation  to  within 
5  percent  at  a  Mach  number  of  0.8  for  thickness  ratios  up  to  0.30.  Although 
the  second  approximation  is  not  the  exact  solution,  it  indicates  that  the 
error  involved  in  using 'the  Prandtl-Glauert  method  to  estimate  the  maximum 
incremental  velocity  for  ellipsoids  having  a  given  thickness  ratio  is 
greatest  for  the  limiting  case  of  the  elliptic  cylinder  (A  =  “)  and 
very  small  for  the  ellipsoid  of  revolution,  which  has  a  very  low  aspect 
ratio.  The  error  may  be  expected  to  be  intermediate  in  ifnagnltnle  for 
intermediate  values  of  the  aspect  ratio  and  to  decrease  with  aspect  ratio. 
The  reduction  of  error  of  the  Prandtl-Glauert  method  with  a  decrease  in 
aspect  ratio  was  to  be  expected,  as  the  Incremental  velocities  are  smaller 
for  ellipsoids  having  low  aspect  ratio. 


RESULTS  MD  DISCUSSION 


Ee suits.  -  Figures  1  to  3  show  the  value  of  the  velocity 


ratio  u  =  — —  at  the  half -chord  line  plotted  against  the  Mach  number 

for  ellipsoids  at  zero  angle  of  attack  for  various  aspect  ratios  and  sec¬ 
tion  thickness  ratios  equal  to  0»10,  0.15^  and  0.20.  In  the  same  figures 
the  sonic  velocity  boundary  having  the  equation 


-  1 


is  plotted  for  air-  The  abscissa  of  the  intersection  of  this  boundary  line 
with  the  curve  of  u  plotted  against  M  for  any  aspect  ratio  is  the 
critical  Mach  number  (within  the  accuracy  of  the  Prandtl-Glauert  method) . 

In  order  to  show  the  effect  of  compressibility  more  directly,  the 
u(m) 

ratio  of  maximum  incremental  velocity  for  compressible  flow  to 


the  maximum  Incremental  velocity  for  incompressible  flow  for  the  same 
free-stream  velocity  is  plotted  against  the  Mach  mmober  in  figures  k 
to  6  for  the  same  aspect  ratios  and  thickness  ratios.  Similar  curves 
for  the  ellipsoid  of  revolution,  which  is  a  special  case  of  the  ellipsoid 
having  three  unequal  axes,  are  plotted  for  the  same  thickness  ratios 
in  figures  1  to  6.  Figure  7  presents  curves  of  critical  Mach  number 
against  aspect  ratio  for  thickness  ratios  of  0-10,  0.15,  and  0.20. 


Three-dimensional  relief.  -  It  may  be  seen  from  figures  1  to  3  that 
the  three-dimensional  relief,  that  is,  the  difference  between  the  velocity 
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on  the  ellipsoid  and  the  velocity  on  the  corresponding  ellipsoid  of  infinite 
aspect  ratio  (elliptic  cylinder),  increases  vith  a  decrease  in  the  aspect 
ratio.  This  increase  has  two  causes: 

(1)  For  a  flow  at  omflU  values  of  M  (incoanpressihle  flow),  the 
relief  effect  increases  with  a  decrease  in  the  aspect  ratio. 

(2)  For  larger  values  of  M  (campressihle  flow),  an  additional 
relief  effect  occurs  with  a  decrease  in  the  aspect  ratio  because  of  the 
fact  that  the  compressibility  effect  (increase  of  incremental  velocity 
with  an  increase  in  the  Mach  number)  decreases  with  a  decrease  in  the 
aspect  ratio.  (See  figs.  4  to  6* )  It  may  be  seen  that  this  additional 
three-dimensional  relief  increases  most  rapidly  at  high  Mach  numbers. 

Figures  1  to  6  show  that  the  compressibility  effect  on  the  maximum 
incremental  velocity  is  greatest  for  A  equal  to  infinity  (infinite 
elliptic  cylinder)  and  is  smallest  for  the  ellipsoid  of  revolution.  The 
compressibility  effect  on  the  maximum  incremental  velocity  for  the  elliptic 

cylinder  is  proportional  to  . — 3 — :,  which  is  in  agreement  with  the 

Vl  - 

usual  form  of  the  Prandtl-Glauert  method  in  two  dimensions.  The  com¬ 
pressibility  effect  on  the  maximum  incremental  velocity  for  the  ellipsoid 
of  revolution  is  small  in  comparison  with  that  of  the  elliptic  cylinder* 

In  fact,  as  the  thickness  ratio  of  any  type  of  body  of  revolution 
approaches  zero,  the  compressibility  correction  factor  approaches  vuilty, 
for  in  this  limit  the  incremental  velocity  in  incompressible  flow  is 
proportional  to  the  square  of  the  thickness  ratio,  so  that  the  effect  of 
stretching  the  body  (first  step  of  Prandtl-Glauert  method,  see  appendix) 
is  exactly  compensated  for  by  the  multiplication  of  the  Incremental 

velocities  by  l/p^  (third  step  of  the  Prandtl-Glauert  method).  For 
ellipsoids  of  practical  thickness  ratios,  however,  the  incremental 
velocity  varies  more  slowly  than  the  square  of  the  thickness  ratio. 

The  compressibility  effect  for  the  ellipsoid  of  revolution  (figs.  4  to  6) 
is  thus  appreciable  at  high  Mach  numbers.  For  example,  for  a  thickness 
ratio  of  0-20  and  at  a  Mach  number  of  0.8,  the  compressibility  effect 
amounts  to  about  30  percent  of  the  Incremental  velocity  in  incompressible 
flow. 


The  effect  of  the  thickness  ratio  on  the  three-dimensional  relief 
may  be  seen  by  a  comparison  of  figures  1,  2,  and  3*  From  figure  1  it 
may  be  seen  that,  for  a  thickness  ratio  of  0.10,  at  a  Mach  number  of  0*75^ 

the  maximum  incremental  velocity  for  A  =  2  is  76  percent  of  the  maximum 

incremental  velocity  for  A  =  From  figure  3^  on  the  other  hand,  it 
may  be  seen  that,  for  a  thickness  ratio  of  0.20,  at  a  Mach  number  of  0.75, 

the  maximum  Incremental  velocity  for  A  =  2  is  75  percent  of  the  maximum 

incremental  velocity  for  A  =  00.  Thus,  an  Increase  in  the  thickness  ratio 
causes  only  a  very  small  Increase  in  the  three-dimensional  relief. 
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Critical  Mach  nuiober*  Figures  1,  2,  3^  and  7  ln-d.lcate  that  an  increase 
in  the  critical  Mach  number  of  an  ellipsoid  at  zero  lift  may  he  obtained 
hy  decreasing  the  aspect  ratio.'  For  example,  for  ellipsoids  having  a 
thickness  ratio  of  0*10,  a  decrease  in  the  aspect  ratio  from  »  to  2 
causes  the  critical  Mach  nimiher  to  increase  from  0-827  to  0.857  (a  Mach 
number  increase  of  O-QS)-  For  a  thickness  ratio  of  0.20,  a  decrease  in 
the  aspect  ratio  from  “  to  2  causes  the  critical  Mach  number  to  increase 
from  0-7^1  to  O.783  (about  0-04).  Although  ellipsoids  having  greater 
thickness  ratio  have  lower  critical  Mach  numbers,  a  decrease  in  the 
aspect  ratio  is  slightly  more  effective  in  increasing  the  critical  Mach 
numbers  for  ellipsoids  of  greater  thickness  ratio.  Figure  7  indicates 
that  only  a  large  reduction  in  aspect  ratio  will  cause  a  significant  rise 
of  the  critical  Mach  number* 

Comuarison  with  test  results  on  low  asnect  ratio  wings.  -  Figure  ^ 
of  reference  1  shows  the  minimum  drag  coefficient  for  zero  lift) 

plotted  against  the  Mach  number  for  wings  having  an  NACA  0012  section 
and  various  aspect  ratios.  The  critical  Mach  number  for  any  aspect 
ratio  may  be  estimated  roughly  as  the  Mach  number  for  which  the  drag 
coefficient  first  begins  to  rise.  The  rou^  estimate  of  the  critical 
Mach  numbers  obtainable  by  this  consideration  is  not  sufficiently  accurate 
to  warrant  comparison  of  the  numerical  values  with  the  numerical  values* 
of  the  critical  Mach  number  obtained  in  the  present  paper  for  thin 
ellipsoids.  Comparison  of  the  numerical  results  is,  moreover,  not 
warranted  inasmuch  as  the  wings  of  reference  1  did  not  have  an  elliptic 
section  and  fuz^thermore  had  a  rectangular  plan  form.  A  q_ualitative 
comparison  may  be  made,  however,  between  the  results  of  the  present 
paper  and  those  of  reference  1.  The  increase  in  critical  Mach  number 
with  decrease  in  aspect  ratio  indicated  in  figures  1,  2,  3,  and  7  of 
the  present  paper  is  considered  sufficiently  large  to  explain  the 
corresponding  effect  indicated  in  figure  6,  reference  1. 

It  is  mentioned  in  reference  1  that  the  Mach  number  for  a  significant 
rise  in  the  drag  coefficient  is  approximately  0*1  hi^er  for  an  aspect 
ratio  of  2  than  for  an  infinite  aspect  ratio.  This  value  is  appreciably 
higher  than  the  increase  in  critical  Mach  number  due  to  a  decrease  in  the 
aspect  ratio.  Since,  for  low-aspect -ratio  wings,  the  drag  coefficient 
increases  only  gradually  after  the  critical  Mach  number  is  reached,  the 
critical  Mach  number  for  a  wing  having  low  aspect  ratio  does  not  indicate 
so  critical  a  change  in  the  flow  phenomena  as  the  critical  Mabh  number 
for  a  wing  having  high  aspect  ratio-  It  is  thought  that  the  smaller  rate 
of  increase  of  the  drag  coefficient  for  wings  having  low  aspect  ratio  is 
due  to  the  fact  that,  at  the  critical  Mach  number,  the  rate  of  increase 
with  Mach  number  of  the  incremental  velocity  is  less  than  for  high  aspect 
ratios,  as  may  be  seen  from  figures  1  to  3* 
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CONCLUSIONS 


A  stud^y-  by  the  Prandtl-Glauert  method  of  compressibility  effects 
and  critical  Mach  number  for  ellipsoids  of  various  aspect  ratios  and 
thickness  ratios  Indicated  the  follo’wrlng  conclusions: 

1.  The  flov  about  the  unyawed  ellipsoid  is  analogous  to  that  about 
the  Infinite  unyawed  cylinder  In  that  sonic  velocity  Is  reached  simul¬ 
taneously  along  a  line  that  extends  across  the  entire  span  of  the  body 
and  Is  normal  to  the  stream  direction. 

2.  The  critical  Mach  number  for  a  thin  ellipsoid  may  be  predicted 
with  good  accuracy  by  means  of  the  Prandtl-Glauert  method,  and  the 
accuracy  increases  with  decrease  in  aspect  ratio. 

3.  The  compressibility  effect  on  the  flow  about  an  ellipsoid  decreases 
as  the  aspect  ratio  decreases. 

4.  The  three-dimensional  relief  for  ellipsoids  is  essentially 
independent  of  the  thickness  ratio,  for  thickness  ratios  from  0*10  to  0.20. 

5.  For  ellipsoids  of  thickness  ratio  0.20,  the  critical  Mach  number 
increases  by  about  0.04  when  the  aspect  ratio  is  changed  from  00  to  2j 
for  ellipsoids  of  thickness  ratio  0.10  the  increase  is  0.03. 

6.  The  calculated  increases  in  critical  Mach  number  are  sufficiently 
large  to  explain  the  experimentally  observed  increases  in  the  Mach  number 
at  which  the  drag  first  begins  to  rise. 

7*  The  experimentally  indicated  reduced  rate  of  drag  rise  for  low- 
aspect-ratio  wings  at  zero  lift  as  compared  with  that  for  wings  having 
infinite  aspect  ratio  may  be  explained  qualitatively  on  the  basis  of  the 
results  obtained  for  the  three-dimensional  relief  for  ellipsoids. 


Langley  Aeronautical  Laboratory 
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APPENDIX 

TEE  PEMDTL-GLAIIERT  METHOD  FOE  THREE-DIMENSIONAL  FLOW 


A  derivation  of  the  Prandtl-Glauert  method  for  three -dlmeiislonal  flov. - 
A  ‘brief  derivation  of  a  form  of  the  Prandtl-Glauert  method  correct  for 
three  dimensions  may  "be  given  as  follows:  A  first-order  approximation  to 
the  subsonic  compressible  flow  about  a  thin  body  B,  the  surface  of  which 
has  the  equation 


S(x,  y,  z)  =  0 

may  be  obtained  by  finding  a  solution  of  the  linearized  differential 
eq^uatlon  for  the  potential  cp  of  the  Incremental  velocities, 


P^^XX  +  'Pyy  +  ^Z2  =  0 

where  the  x-axls  is  in  the  stream  direction  and  the  Incremental 

velocities  cp  ^  and  cp  are  small  compared  with  the  stream 

^  y  z 

velocity  IJ.  At  all  points  on  the  surface  of  B,  the  potential  cp  must 
satisfy  the  houndary.  condition 

(u  +  tp^)  S^  +  cPySy  +  =  0  (A2) 

which  states  that  the  flow  is  tangential  to  B.  Since  B  is  assumed 
thin,  Sx  is  small '  compared  with  Sj  and  S^J  conseq^uently  the  second- 

order  teim  cp  s  may  he  neglected,  and  the  houndary  condition  becomes 


In  order  to  solve  the  houndary -value  problem  given  by  equations  (Al) 
and  (A2)  in  terms  of  incompressible  flow,  the  following  transformation  of 
variables  is  used: 


X 


! 


q)' 


(A3) 
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ll+ 


With  this  transformation  equations  (Al)  and  (A2)  become,  respectively. 


<p'  1  (  +  cp' +  cp'  =0 

X  X  ^  yy  ^  zz 

.  (aA) 

US  ,  +  cp '  s  +  cp  '  =  0 

X  ’  '  y  y  ^  z  z 

(A5) 

Equations  (AU)  and  (A5)  are,  respectively,  the  differential  equation  and 
boundary  condition  for  the  potential  cp'  of  the  incremental  velocities 
of  an  incompressible  flow  with  free^stream  velocity  U,  in  the  x',  y,  z 
space,  about  a  thin  body  B',  the  surface  of  which  has  the  equation 


S(Px',  y,  z)  =  0 

The  incremental  velocities  in  the  compressible  flow  are  thus 
given  by 

u  =  cp^  =  ^  cp'x'  =  ^ 

1  .  1  I 

where  u,  v,  and  w  and  v’,  and  w*  are  the  incremental 

velocities  at  corresponding  points  in  the  compressible  flow  about  B 
and  the  incompressible  flow  about  B',  respectively. 

The  foregoing  analysis  establishes  the  Prandtl-Glauert  method 
for  three-dimensional  flow  in  the  following  form:  The  incremental 
velocities  at  a  point  P  of  a  three-dimensional  compressible  flow 
field  about  a  thin  body  B  may  be  obtained  in  three  steps: 

(1)  The  x-coordinates  of  all  points  of  B  are  Increased  by  the 
factor  1/3,  where 

P  =  ^1  -  m2 

and  where  the  x-axis  is  in  the  stream  direction.  This  transformation 
takes  B  into  a  stretched  body  B*. 

(2)  The  incremental  velocities  u',  v',  w*,  in  the  direction  of 
the  X-;  y-,  and  z-axes,  respectively,  at  the  point  P'  in  the  flow 
field  of  B'  corresponding  to  the  point  P  in  the  flow  field  of  B 
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are  calculated  as  though  B '  were  in  an  incompressihle  flow  having  the 
same  free-stream  velocity  as  the  original  compressihle  flow. 

(3)  The  values  u,  v,  and  w  of  the  Incremental  velocities  at  the 
point  P  in  the  compressihle  flow  field  of  the  original  unstretched 
body  are  then  found  hy  the  following  equations: 


u 


V 


1. 

32 

P 


u 


I 


V 


t 


Thus  far  It  has  heen  shovn  that^  through  the  transformation  given 
in  this  paper,  a  compressible  flov  that  satisfies  the  boundary  conditions 
for  the  body  B  is  transformed  into  an  incompressible  flow  satisfying 
the  boundary  conditions  for  the  stretched  body  B*.  It  can  be  shown 
further  that  the  stream  lines  of  the  compressible  flow  about  B  are 
transformed  into  stream  lines  of  the  Incompressible  flow  about  B*. 
Because  of  this  fact,  the  method  has  been  referred  to  in  the  literature 
(for  example,  reference  4)  as  the  "streamline -analogy method. 


The  proof  is  obtained  simply  by  applying  the  transformations  (A3) 
and  (a6)  to  the  ec[uatlons  for  the  streamlines  of  the  incompressible 
flow  about  the  stretched  body  B* 


dx '  dy  ‘  dz  ' 
U  v^  w’ 


(AT) 


The  application  of  the  transformations  resxiLts  in  the  equation  of 
the  streamlines  for  the  compressible  flow  about  the  body  B 


U  V  w 


(a8) 


Failure  for  three-dimensional  flow  -problems  of  the  commonly  stated 
foimis  of  the  Prandtl-Qlauert  method.  -  According  to  ^he  form  of  the 
Prandtl-Glauert  method  given  by  Prandtl  (reference  5)  and  Yon  Karm^ 
(reference  6),  the  incremental  velocities  for  a  compressible  flow  about 
a  thin  body  B  are  the  same  as  the  incremental  velocities  of  corresponding 
points  for  Incompresalble  flow  having  the  same  free-stream  velocity  about 
a  body  obtained  by  expanding  B  in  the  directions  normal  to  the  free- 
stream  direction  by  the  factor  l/p.  That  is,  for  bodies  of  revolution, 
or  two-dimensional  bodies. 
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According  to  Gothert ' s  method,  hovever , 

^  (e,  M)  =  ^  Ji  Oe,  0)  (A9) 

u  u 

Thus,  Prandtl's  and  Von  Karmm's  method  is  valid  only  if 

that  is.  If  and  only  If  the  Incremental  velocity  for  incompressihle 
flow  about  the  bodies  under  consideration  is  proportional  to  the  thickness 
ratio.  This  relation  is  approximately  valid  for  thin  two-dimensional 
bodies,  so  that  the  method  of  Prandtl  and  Von  Karmen  may  be  expected 
to  be  vsdid  for  two-dimensional  flows.  The  relation  is  not  true  in 
general  for  three-dimensional  bodies;  for  example,  for  a  very  thin 
body  of  revolution  the  Incremental  velocity  is  more  nearly  proportional 
to  the  square  of  the  thickness  ratio  than  to  the  first  power. 

Von  approaches  the  problem  by  making  the  transformation 

y'  =  3y 

z  '  =  Pz 
cp'  =  cp 

Under  this  transformation  the  linearized  equation  of  compressible 
flow  goes  into  Laplace's  equation;  however,  the  transformed  boundary 
condition  is  not  satisfied  on  the  siirface  of  the  transformed  (contracted) 
body  but  on  the  surface  of  an  expanded  body.  Thus,  the  boundary  condi¬ 
tion  is  not  satisfied  on  the  boundary  but  at  points  near  the  boundary. 

This  procedure  is  applicable  to  two-dimensional  problems  (as,  for  example, 
in  the  thln-wing  theory,  reference  14) ,  because  the  velocity  increments 
Induced  by  the  equivalent  line  distribution  of  singularities  vary  only 
slowly  in  the  neighborhood  of  the  line  of  singularities.  For  a  body 
of  revolution,  however,  the  velocity  increments  Induced  by  a  line  of 
singularities  go  to  infinity  at  the  line  of  singularities;  for  such 
bodies,  accordingly,  the  location  of  the  point  at  which  the  boundaiy 
condition  is  satisfied  is  important. 

According  to  Goldstein  and  Young  (reference  7)^  "in  compressible 
flow  the  pressure  increase  at  any  point  of  the  body  is  l/p  times  the 
pressure  Increase  in  Incompressible  flow  at  the  same  point.  That  is. 
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M)  =71^.0) 


Comparison  of  this  relation  with  equation  (A9)  shows  that  the  Goldstein- 
Young  method  Is  also  valid  for  two-dimensional  problems  hut  gives  an 
Incorrect  result  for  three-dimensional  problems- 
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Hgure  3-  Variation  w/th  Mach  number  of  maxi  mam  incremental  vekxitu  for 
ellipsoids  hamg  various  aspect  ratios.  Thickness  ratio.  0.20. 


rigure4-  Variation  of  cgmoressibiUty  effect  w/th  Mach  number  for 
ellipsoids  ha\/ing  various  aspect  ratios.  Thickness  ratio ^QIO. 


Figure  6.-  Variation  of  compressibility  yffect  vyift)  Fach  number  for  ellipsoids 
having  various  aspect  ratios^  Thickness  ratio,  0.2.0. 


/V  ’^dgujnu  / 


Figure  7.  -  Variat/on  of  gritica/  Mach  number  with  aspect  ratio 
el/ipsoips  hay/ng  thickness  ratios  of  OJO^  O./b  and  0.£ 


